Porcine deltacoronavirus accessory protein NS7a possesses the functional characteristics of a viroporin.

Viroporins are virus-encoded proteins that mediate ion channel (IC) activity, playing critical roles in virus entry, replication, pathogenesis, and immune evasion. Previous studies have shown that some coronavirus accessory proteins have viroporin-like activity. Porcine deltacoronavirus (PDCoV) is an emerging enteropathogenic coronavirus that encodes three accessory proteins, NS6, NS7, and NS7a. However, whether any of the PDCoV accessory proteins possess viroporin-like activity, and if so which, remains unknown. In this study, we analyzed the biochemical properties of the three PDCoV-encoded accessory proteins and found that NS7a could enhance the membrane permeability of both mammalian cells and Escherichia coli cells. Indirect immunofluorescence assay and co-immunoprecipitation assay results further indicated that NS7a is an integral membrane protein and can form homo-oligomers. A bioinformation analysis revealed that a putative viroporin domain (VPD) is located within amino acids 69-88 (aa69-88) of NS7a. Experiments with truncated mutants and alanine scanning mutagenesis additionally demonstrated that the amino acid residues 69FLR71 of NS7a are essential for its viroporin-like activity. Together, our findings are the first to demonstrate that PDCoV NS7a possesses viroporin-like activity and identify its key amino acid residues associated with viroporin-like activity.

In Vitro Selection of Remdesivir-Resistant SARS-CoV-2 Demonstrates High Barrier to Resistance.

In vitro selection of remdesivir-resistant severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) revealed the emergence of a V166L substitution, located outside of the polymerase active site of the Nsp12 protein, after 9 passages of a single lineage. V166L remained the only Nsp12 substitution after 17 passages (10 µM remdesivir), conferring a 2.3-fold increase in 50% effective concentration (EC50). When V166L was introduced into a recombinant SARS-CoV-2 virus, a 1.5-fold increase in EC50 was observed, indicating a high in vitro barrier to remdesivir resistance.

Mutational Effect of Some Major COVID-19 Variants on Binding of the S Protein to ACE2.

COVID-19 is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which has many variants that accelerated the spread of the virus. In this study, we investigated the quantitative effect of some major mutants of the spike protein of SARS-CoV-2 binding to the human angiotensin-converting enzyme 2 (ACE2). These mutations are directly related to the Variant of Concern (VOC) including Alpha, Beta, Gamma, Delta and Omicron. Our calculations show that five major mutations (N501Y, E484K, L452R, T478K and K417N), first reported in Alpha, Beta, Gamma and Delta variants, all increase the binding of the S protein to ACE2 (except K417N), consistent with the experimental findings. We also studied an additional eight mutations of the Omicron variant that are located on the interface of the receptor binding domain (RDB) and have not been reported in other VOCs. Our study showed that most of these mutations (except Y505H and G446S) enhance the binding of the S protein to ACE2. The computational predictions helped explain why the Omicron variant quickly became dominant worldwide. Finally, comparison of several different computational methods for binding free energy calculation of these mutants was made. The alanine scanning method used in the current calculation helped to elucidate the residue-specific interactions responsible for the enhanced binding affinities of the mutants. The results show that the ASGB (alanine scanning with generalized Born) method is an efficient and reliable method for these binding free energy calculations due to mutations.

Comparative assessment of favipiravir and remdesivir against human coronavirus NL63 in molecular docking and cell culture models.

Human coronavirus NL63 (HCoV-NL63) mainly affects young children and immunocompromised patients, causing morbidity and mortality in a subset of patients. Since no specific treatment is available, this study aims to explore the anti-SARS-CoV-2 agents including favipiravir and remdesivir for treating HCoV-NL63 infection. We first successfully modelled the 3D structure of HCoV-NL63 RNA-dependent RNA polymerase (RdRp) based on the experimentally solved SARS-CoV-2 RdRp structure. Molecular docking indicated that favipiravir has similar binding affinities to SARS-CoV-2 and HCoV-NL63 RdRp with LibDock scores of 75 and 74, respectively. The LibDock scores of remdesivir to SARS-CoV-2 and HCoV-NL63 were 135 and 151, suggesting that remdesivir may have a higher affinity to HCoV-NL63 compared to SARS-CoV-2 RdRp. In cell culture models infected with HCoV-NL63, both favipiravir and remdesivir significantly inhibited viral replication and production of infectious viruses. Overall, remdesivir compared to favipiravir is more potent in inhibiting HCoV-NL63 in cell culture. Importantly, there is no evidence of resistance development upon long-term exposure to remdesivir. Furthermore, combining favipiravir or remdesivir with the clinically used antiviral cytokine interferon-alpha resulted in synergistic effects. These findings provided a proof-of-concept that anti-SARS-CoV-2 drugs, in particular remdesivir, have the potential to be repurposed for treating HCoV-NL63 infection.

Human carboxylesterase 1A plays a predominant role in the hydrolytic activation of remdesivir in humans.

Remdesivir, an intravenous nucleotide prodrug, has been approved for treating COVID-19 in hospitalized adults and pediatric patients. Upon administration, remdesivir can be readily hydrolyzed to form its active form GS-441524, while the cleavage of the carboxylic ester into GS-704277 is the first step for remdesivir activation. This study aims to assign the key enzymes responsible for remdesivir hydrolysis in humans, as well as to investigate the kinetics of remdesivir hydrolysis in various enzyme sources. The results showed that remdesivir could be hydrolyzed to form GS-704277 in human plasma and the microsomes from human liver (HLMs), lung (HLuMs) and kidney (HKMs), while the hydrolytic rate of remdesivir in HLMs was the fastest. Chemical inhibition and reaction phenotyping assays suggested that human carboxylesterase 1 (hCES1A) played a predominant role in remdesivir hydrolysis, while cathepsin A (CTSA), acetylcholinesterase (AchE) and butyrylcholinesterase (BchE) contributed to a lesser extent. Enzymatic kinetic analyses demonstrated that remdesivir hydrolysis in hCES1A (SHUTCM) and HLMs showed similar kinetic plots and much closed Km values to each other. Meanwhile, GS-704277 formation rates were strongly correlated with the CES1A activities in HLM samples from different individual donors. Further investigation revealed that simvastatin (a therapeutic agent for adjuvant treating COVID-19) strongly inhibited remdesivir hydrolysis in both recombinant hCES1A and HLMs. Collectively, our findings reveal that hCES1A plays a predominant role in remdesivir hydrolysis in humans, which are very helpful for predicting inter-individual variability in response to remdesivir and for guiding the rational use of this anti-COVID-19 agent in clinical settings.

Remdesivir: Quo vadis?

Remdesivir (GS-5734, Veklury®) has remained the only antiviral drug formally approved by the US FDA for the treatment of Covid-19 (SARS-CoV-2 infection). Its key structural features are the fact that it is a C-nucleoside (adenosine) analogue, contains a 1'-cyano function, and could be considered as a ProTide based on the presence of a phosphoramidate group. Its antiviral spectrum and activity in animal models have been well established and so has been its molecular mode of action as a delayed chain terminator of the viral RdRp (RNA-dependent RNA polymerase). Its clinical efficacy has been evaluated, but needs to be optimized with regard to timing, dosage and duration of treatment, and route of administration. Safety, toxicity and pharmacokinetics need to be further addressed, and so are its potential combinations with other drugs such as corticosteroids (i.e. dexamethasone) and ribavirin.

Whole Body PBPK Modeling of Remdesivir and Its Metabolites to Aid in Estimating Active Metabolite Exposure in the Lung and Liver in Patients With Organ Dysfunction.

Remdesivir (RDV) is the first drug approved by the US Food and Drug Administration (FDA) for the treatment of coronavirus disease 2019 (COVID-19) in certain patients requiring hospitalization. As a nucleoside analogue prodrug, RDV undergoes intracellular multistep activation to form its pharmacologically active species, GS-443902, which is not detectable in the plasma. A question arises that whether the observed plasma exposure of RDV and its metabolites would correlate with or be informative about the exposure of GS-443902 in tissues. A whole body physiologically-based pharmacokinetic (PBPK) modeling and simulation approach was utilized to elucidate the disposition mechanism of RDV and its metabolites in the lungs and liver and explore the relationship between plasma and tissue pharmacokinetics (PK) of RDV and its metabolites in healthy subjects. In addition, the potential alteration of plasma and tissue PK of RDV and its metabolites in patients with organ dysfunction was explored. Our simulation results indicated that intracellular exposure of GS-443902 was decreased in the liver and increased in the lungs in subjects with hepatic impairment relative to the subjects with normal liver function. In subjects with severe renal impairment, the exposure of GS-443902 in the liver was slightly increased, whereas the lung exposure of GS-443902 was not impacted. These predictions along with the organ impairment study results may be used to support decision making regarding the RDV dosage adjustment in these patient subgroups. The modeling exercise illustrated the potential of whole body PBPK modeling to aid in decision making for nucleotide analogue prodrugs, particularly when the active metabolite exposure in the target tissues is not available.

Effect of drug metabolism in the treatment of SARS-CoV-2 from an entirely computational perspective.

Understanding the effects of metabolism on the rational design of novel and more effective drugs is still a considerable challenge. To the best of our knowledge, there are no entirely computational strategies that make it possible to predict these effects. From this perspective, the development of such methodologies could contribute to significantly reduce the side effects of medicines, leading to the emergence of more effective and safer drugs. Thereby, in this study, our strategy is based on simulating the electron ionization mass spectrometry (EI-MS) fragmentation of the drug molecules and combined with molecular docking and ADMET models in two different situations. In the first model, the drug is docked without considering the possible metabolic effects. In the second model, each of the intermediates from the EI-MS results is docked, and metabolism occurs before the drug accesses the biological target. As a proof of concept, in this work, we investigate the main antiviral drugs used in clinical research to treat COVID-19. As a result, our strategy made it possible to assess the biological activity and toxicity of all potential by-products. We believed that our findings provide new chemical insights that can benefit the rational development of novel drugs in the future.

Indole Chloropyridinyl Ester-Derived SARS-CoV-2 3CLpro Inhibitors: Enzyme Inhibition, Antiviral Efficacy, Structure-Activity Relationship, and X-ray Structural Studies.

Here, we report the synthesis, structure-activity relationship studies, enzyme inhibition, antiviral activity, and X-ray crystallographic studies of 5-chloropyridinyl indole carboxylate derivatives as a potent class of SARS-CoV-2 chymotrypsin-like protease inhibitors. Compound 1 exhibited a SARS-CoV-2 3CLpro inhibitory IC50 value of 250 nM and an antiviral EC50 value of 2.8 µM in VeroE6 cells. Remdesivir, an RNA-dependent RNA polymerase inhibitor, showed an antiviral EC50 value of 1.2 µM in the same assay. Compound 1 showed comparable antiviral activity with remdesivir in immunocytochemistry assays. Compound 7d with an N-allyl derivative showed the most potent enzyme inhibitory IC50 value of 73 nM. To obtain molecular insight into the binding properties of these molecules, X-ray crystal structures of compounds 2, 7b, and 9d-bound to SARS-CoV 3CLpro were determined, and their binding properties were compared.

Remdesivir and EIDD-1931 Interact with Human Equilibrative Nucleoside Transporters 1 and 2: Implications for Reaching SARS-CoV-2 Viral Sanctuary Sites.

Equilibrative nucleoside transporters (ENTs) are present at the blood-testis barrier (BTB), where they can facilitate antiviral drug disposition to eliminate a sanctuary site for viruses detectable in semen. The purpose of this study was to investigate ENT-drug interactions with three nucleoside analogs, remdesivir, molnupiravir, and molnupiravir's active metabolite, ß-d-N4-hydroxycytidine (EIDD-1931), and four non-nucleoside molecules repurposed as antivirals for coronavirus disease 2019 (COVID-19). The study used three-dimensional pharmacophores for ENT1 and ENT2 substrates and inhibitors and Bayesian machine learning models to identify potential interactions with these transporters. In vitro transport experiments demonstrated that remdesivir was the most potent inhibitor of ENT-mediated [3H]uridine uptake (ENT1 IC50: 39 µM; ENT2 IC50: 77 µM), followed by EIDD-1931 (ENT1 IC50: 259 µM; ENT2 IC50: 467 µM), whereas molnupiravir was a modest inhibitor (ENT1 IC50: 701 µM; ENT2 IC50: 851 µM). Other proposed antivirals failed to inhibit ENT-mediated [3H]uridine uptake below 1 mM. Remdesivir accumulation decreased in the presence of 6-S-[(4-nitrophenyl)methyl]-6-thioinosine (NBMPR) by 30% in ENT1 cells (P = 0.0248) and 27% in ENT2 cells (P = 0.0054). EIDD-1931 accumulation decreased in the presence of NBMPR by 77% in ENT1 cells (P = 0.0463) and by 64% in ENT2 cells (P = 0.0132), which supported computational predictions that both are ENT substrates that may be important for efficacy against COVID-19. NBMPR failed to decrease molnupiravir uptake, suggesting that ENT interaction is likely inhibitory. Our combined computational and in vitro data can be used to identify additional ENT-drug interactions to improve our understanding of drugs that can circumvent the BTB. SIGNIFICANCE STATEMENT: This study identified remdesivir and EIDD-1931 as substrates of equilibrative nucleoside transporters 1 and 2. This provides a potential mechanism for uptake of these drugs into cells and may be important for antiviral potential in the testes and other tissues expressing these transporters.

Interactions of anti-COVID-19 drug candidates with hepatic transporters may cause liver toxicity and affect pharmacokinetics.

Transporters in the human liver play a major role in the clearance of endo- and xenobiotics. Apical (canalicular) transporters extrude compounds to the bile, while basolateral hepatocyte transporters promote the uptake of, or expel, various compounds from/into the venous blood stream. In the present work we have examined the in vitro interactions of some key repurposed drugs advocated to treat COVID-19 (lopinavir, ritonavir, ivermectin, remdesivir and favipiravir), with the key drug transporters of hepatocytes. These transporters included ABCB11/BSEP, ABCC2/MRP2, and SLC47A1/MATE1 in the canalicular membrane, as well as ABCC3/MRP3, ABCC4/MRP4, SLC22A1/OCT1, SLCO1B1/OATP1B1, SLCO1B3/OATP1B3, and SLC10A1/NTCP, residing in the basolateral membrane. Lopinavir and ritonavir in low micromolar concentrations inhibited BSEP and MATE1 exporters, as well as OATP1B1/1B3 uptake transporters. Ritonavir had a similar inhibitory pattern, also inhibiting OCT1. Remdesivir strongly inhibited MRP4, OATP1B1/1B3, MATE1 and OCT1. Favipiravir had no significant effect on any of these transporters. Since both general drug metabolism and drug-induced liver toxicity are strongly dependent on the functioning of these transporters, the various interactions reported here may have important clinical relevance in the drug treatment of this viral disease and the existing co-morbidities.

Mechanism of Inhibition of the Reproduction of SARS-CoV-2 and Ebola Viruses by Remdesivir.

Remdesivir is an antiviral drug initially designed against the Ebola virus. The results obtained with it both in biochemical studies in vitro and in cell line assays in vivo were very promising, but it proved to be ineffective in clinical trials. Remdesivir exhibited far better efficacy when repurposed against SARS-CoV-2. The chemistry that accounts for this difference is the subject of this study. Here, we examine the hypothesis that remdesivir monophosphate (RMP)-containing RNA functions as a template at the polymerase site for the second run of RNA synthesis, and as mRNA at the decoding center for protein synthesis. Our hypothesis is supported by the observation that RMP can be incorporated into RNA by the RNA-dependent RNA polymerases (RdRps) of both viruses, although some of the incorporated RMPs are subsequently removed by exoribonucleases. Furthermore, our hypothesis is consistent with the fact that RdRp of SARS-CoV-2 selects RMP for incorporation over AMP by 3-fold in vitro, and that RMP-added RNA can be rapidly extended, even though primer extension is often paused when the added RMP is translocated at the i + 3 position (with i the nascent base pair at an initial insertion site of RMP) or when the concentrations of the subsequent nucleoside triphosphates (NTPs) are below their physiological concentrations. These observations have led to the hypothesis that remdesivir might be a delayed chain terminator. However, that hypothesis is challenged under physiological concentrations of NTPs by the observation that approximately three-quarters of RNA products efficiently overrun the pause.

The Enzyme-Free Release of Nucleotides from Phosphoramidates Depends Strongly on the Amino Acid.

Phosphoramidates composed of an amino acid and a nucleotide analogue are critical metabolites of prodrugs, such as remdesivir. Hydrolysis of the phosphoramidate liberates the nucleotide, which can then be phosphorylated to become the pharmacologically active triphosphate. Enzymatic hydrolysis has been demonstrated, but a spontaneous chemical process may also occur. We measured the rate of enzyme-free hydrolysis for 17 phosphoramidates of ribonucleotides with amino acids or related compounds at pH 7.5. Phosphoramidates of proline hydrolyzed fast, with a half-life time as short as 2.4 h for Pro-AMP in ethylimidazole-containing buffer at 37 °C; 45-fold faster than Ala-AMP and 120-fold faster than Phe-AMP. Crystal structures of Gly-AMP, Pro-AMP, ßPro-AMP and Phe-AMP bound to RNase A as crystallization chaperone showed how well the carboxylate is poised to attack the phosphoramidate, helping to explain this reactivity. Our results are significant for the design of new antiviral prodrugs.

A fast protein binding site comparison algorithm for proteome-wide protein function prediction and drug repurposing.

The expansion of three-dimensional protein structures and enhanced computing power have significantly facilitated our understanding of protein sequence/structure/function relationships. A challenge in structural genomics is to predict the function of uncharacterized proteins. Protein function deconvolution based on global sequence or structural homology is impracticable when a protein relates to no other proteins with known function, and in such cases, functional relationships can be established by detecting their local ligand binding site similarity. Here, we introduce a sequence order-independent comparison algorithm, PocketShape, for structural proteome-wide exploration of protein functional site by fully considering the geometry of the backbones, orientation of the sidechains, and physiochemical properties of the pocket-lining residues. PocketShape is efficient in distinguishing similar from dissimilar ligand binding site pairs by retrieving 99.3% of the similar pairs while rejecting 100% of the dissimilar pairs on a dataset containing 1538 binding site pairs. This method successfully classifies 83 enzyme structures with diverse functions into 12 clusters, which is highly in accordance with the actual structural classification of proteins classification. PocketShape also achieves superior performances than other methods in protein profiling based on experimental data. Potential new applications for representative SARS-CoV-2 drugs Remdesivir and 11a are predicted. The high accuracy and time-efficient characteristics of PocketShape will undoubtedly make it a promising complementary tool for proteome-wide protein function inference and drug repurposing study.

Mapping the effect of drugs on ACE2 as a novel target site for COVID-19 therapy.

Angiotensin converting enzyme 2 (ACE2) has potentially conflicting roles in health and disease. COVID-19 coronavirus binds to human cells via ACE2 receptor, which is expressed on almost all body organs. Boosting the ACE2 receptor levels on heart and lung cells may provide more cellular enter to virus thereby worsening the infection. Therefore, among the drug targets, ACE2 is suggested as a vital target of COVID-19 therapy. This hypothesis is based on the protective role of the drugs acting on ACE2. Therefore, this review discusses the impact and challenges of using ACE2 as a target in the current therapy of COVID-19.

Role of OATP4C1 in Renal Handling of Remdesivir and its Nucleoside Analog GS-441524: The First Approved Drug for Patients with COVID-19.

PURPOSE: Remdesivir and its active metabolite are predominantly eliminated via renal route; however, information regarding renal uptake transporters is limited. In the present study, the interaction of remdesivir and its nucleoside analog GS-441524 with OATP4C1 was evaluated to provide the detailed information about its renal handling. METHODS: We used HK-2 cells, a proximal tubular cell line derived from normal kidney, to confirm the transport of remdesivir and GS-441524. To assess the involvement of OATP4C1 in handling remdesivir and GS-441524, the uptake study of remdesivir and GS-441524 was performed by using OATP4C1-overexpressing Madin-Darby canine kidney II (MDCKII) cells. Moreover, we also evaluated the IC50 and Ki value of remdesivir. RESULTS: The time-dependent remdesivir uptake in HK-2 cells was observed. The results of inhibition study using OATs and OCT2 inhibitors and OATP4C1 knockdown suggested the involvement of renal drug transporter OATP4C1. Remdesivir was taken up by OATP4C1/MDCKII cells. OATP4C1-mediated uptake of remdesivir increased linearly up to 10 min and reached a steady state at 30 min. Remdesivir inhibited OATP4C1-mediated transport in a concentration-dependent manner with the IC50 and apparent Ki values of 42 ± 7.8 µM and 37 ± 6.9 µM, respectively. CONCLUSIONS: We have provided novel information about renal handling of remdesivir. Furthermore, we evaluated the potential drug interaction via OATP4C1 by calculating the Ki value of remdesivir. OATP4C1 may play a pivotal role in remdesivir therapy for COVID-19, particularly in patients with kidney injury.

Synthesis, molecular docking, and in silico ADMET studies of 4-benzyl-1-(2,4,6-trimethyl-benzyl)-piperidine: Potential Inhibitor of SARS-CoV2.

Nowadays, over 200 countries face a wellbeing emergency because of epidemiological disease COVID-19 caused by the SARS-CoV-2 virus. It will cause a very high effect on the world's economy and the worldwide health sector. The present work is an investigation of the newly synthesized 4-benzyl-1-(2,4,6-trimethyl-benzyl)-piperidine (M1BZP) molecule's inhibitory potential against important protein targets of SARS-CoV-2 using computational approaches. M1BZP crystallizes in monoclinic type with P1211 space group. For the title compound M1BZP, spectroscopic characterization like 1H NMR, 13C NMR, FTIR, were carried out. The geometry of the compound had been optimized by the DFT method and its results were compared with the X-ray diffraction data. The calculated energies for the Highest Occupied Molecular Orbital (HOMO) and the Lowest Unoccupied Molecular Orbital (LUMO) showed the stability and reactivity of the title compound. Intermolecular interactions in the crystal network were determined using Hirshfeld surface analyses. The molecular electrostatic potential (MEP) picture was drawn using the same level of theory to visualize the chemical reactivity and charge distribution on the molecule. Molecular docking study performed for the synthesized compound revealed an efficient interaction with the COVID-19 protease and resulted in good activities. We hope the present study would help workers in the field to develop potential vaccines and therapeutics against the novel coronavirus. Virtual ADME studies were carried out as well and a relationship between biological, electronic, and physicochemical qualifications of the target compound was determined. Toxicity prediction by computational technique for the title compound was also carried out.

Computational hunting of natural active compounds as an alternative for Remdesivir to target RNA-dependent polymerase.

The hunt for potential lead/drug molecules from different resources, especially from natural resources, for possible treatment of COVID-19 is ongoing. Several compounds have already been identified, but only a few are good enough to show potential against the virus. Among the identified druggable target proteins of SARS-CoV-2, this study focuses on non-structural RNA-dependent RNA polymerase protein (RdRp), a well-known enzyme for both viral genome replication and viral mRNA synthesis, and is therefore considered to be the primary target. In this study, the virtual screening followed by an in-depth docking study of the Compounds Library found that natural compound Cyclocurcumin and Silybin B have strong interaction with RdRp and much better than the remdesivir with free binding energy and inhibition constant value as êzŒ-6.29 kcal/mol and 58.39 µMêzŒ, and êzŒ-7.93kcal/mol and 45.3 µMêzŒ, respectively. The finding indicated that the selected hits (Cyclocurcumin and Silybin B) could act as non-nucleotide anti-polymerase agents, and can be further optimized as a potential inhibitor of RdRp by benchwork experiments.

1'-Ribose cyano substitution allows Remdesivir to effectively inhibit nucleotide addition and proofreading during SARS-CoV-2 viral RNA replication.

COVID-19 has recently caused a global health crisis and an effective interventional therapy is urgently needed. Remdesivir is one effective inhibitor for SARS-CoV-2 viral RNA replication. It supersedes other NTP analogues because it not only terminates the polymerization activity of RNA-dependent RNA polymerase (RdRp), but also inhibits the proofreading activity of intrinsic exoribonuclease (ExoN). Even though the static structure of Remdesivir binding to RdRp has been solved and biochemical experiments have suggested it to be a "delayed chain terminator", the underlying molecular mechanisms is not fully understood. Here, we performed all-atom molecular dynamics (MD) simulations with an accumulated simulation time of 24 microseconds to elucidate the inhibitory mechanism of Remdesivir on nucleotide addition and proofreading. We found that when Remdesivir locates at an upstream site in RdRp, the 1'-cyano group experiences electrostatic interactions with a salt bridge (Asp865-Lys593), which subsequently halts translocation. Our findings can supplement the current understanding of the delayed chain termination exerted by Remdesivir and provide an alternative molecular explanation about Remdesivir's inhibitory mechanism. Such inhibition also reduces the likelihood of Remdesivir to be cleaved by ExoN acting on 3'-terminal nucleotides. Furthermore, our study also suggests that Remdesivir's 1'-cyano group can disrupt the cleavage site of ExoN via steric interactions, leading to a further reduction in the cleavage efficiency. Our work provides plausible and novel mechanisms at the molecular level of how Remdesivir inhibits viral RNA replication, and our findings may guide rational design for new treatments of COVID-19 targeting viral replication.